Entanglement sum rules in exactly solvable models 
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We compute the entanglement entropy of a wide class of exactly solvable models which may be characterized 
as describing matter coupled to gauge fields. Our principle result is an entanglement sum rule which states 
that entropy of the full system is the sum of the entropies of the two components. In the context of the exactly 
solvable models we consider, this result applies to the full entropy, but more generally it is a statement about the 
additivity of universal terms in the entropy. We also prove that the Renyi entropy is exactly additive and hence 
that the entanglement spectrum factorizes. Our proof simultaneously extends and simplifies previous arguments, 
with extensions including new models at zero temperature as well as the ability to treat finite temperature 
crossovers. We emphasize that while the additivity is an exact statement, each term in the sum may still be 
difficult to compute. Our results apply to a wide variety of phases including Fermi liquids, spin liquids, and 
some non-Fermi liquid metals. 



I. INTRODUCTION 

A recent exchange of ideas between quantum many-body 
physics and quantum information science has led to an in- 
creased appreciation for the fundamental role of entanglement 
in quantum matter. In particular, long range entanglement 
underlies many of the more interesting states of matter now 
known experimentally, including Fermi liquids uJ4 3[|, q uan- 
tum critical points Jj-[al, and topological phases Il9l- [l3ll . In 
some cases, entanglement is essentially the only completely 
general probe of such states Bl4 Il5ll and has led to a clear 
identification of a topological phase. Entanglement considera- 
tions have also led to a variety of other results, including a new 
class of variational statej K|- T^] and a classification of phases 
in one dimension IU9l 12011 . The concept of entanglement en- 
tropy has played a crucial role in these recent developments, 
so we first remind the reader about entanglement entropy. 

We consider a large quantum system divided into two parts, 
A and B, such that the whole system is in a pure state \4>ab) 
e.g. the ground state of a local Hamiltonian. The entangle- 
ment entropy of A is defined as the von Neumann entropy of 
the reduced density matrix of A: S(A) = — tr a(pa In (pa))- 
Of course this definition makes sense in general, but only 
when the total system is pure does the entanglement entropy 
truly measure entanglement between A and B. Typically A 
and B are spatial regions, but other kinds of entanglement cuts 
have been considered. It is also useful to consider the Renyi 
(entanglement) entropy, defined as 



S n (A) 
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ln(tr(pl)), 



(1.1) 



as a generalization of entanglement entropy. Knowledge of 
the Renyi entropy for all n is equivalent to knowing the full 
spectrum of p^. 

The basic fact about entanglement entropy in local ground 
states is the area law l2ll l22ll . This law states that the entan- 
glement entropy typically scales like S(A) ~ L d_1 where L 
is the linear size of A and d is the spatial dimension. How- 
ever, it must be immediately emphasized that the area law is 
not completely universal as exceptions known in a variety of 
gapless systems. Conversely, although it is not proven, it is be- 
lieved that the area law holds for all gapped phases of matter 



(see Ref. lT23ll for an renormalization group argument and Ref. 
H24I1 for a partial result). The exceptions include conformal 
field theories in one dimension where S ~ cln (L) (c is the 
central charge) [4, 5] and Fermi liquids in d > 1 dimensions 
where S ~ L d_1 In (L) (as well as other systems with a Fermi 
surface) HrH, [25, 26]. Most other known gapped and gapless 
phases in d > 1 satisfy the area law. This central feature of 
entanglement entropy has already led to a new class of varia- 
tional states designed to encode the area law and which extend 
the powerful numerical method known as DMRG i27l l28tl . In 
addition, entanglement considerations have led to a fingerprint 
for topological phases and other strongly correlated systems 
as well as classification schemes for gapped quantum matter. 

However, despite these many advances, entanglement re- 
mains fundamentally poorly understood, especially in the con- 
text of gapless systems. Part of this state of affairs stems from 
the great difficulty, both theoretical and experimental, encoun- 
tered when trying to compute or measure the entropy. Fortu- 
nately, there are a number of promising new directions that 
have recently considerably expanded our ability to compute 
entropies and given us new numerical tests. Here we focus 
on what may be called entanglement sum rules. An entan- 
glement sum rule provides a way to compute the entangle- 
ment properties of an interesting phase of matter by dividing 
it into more elementary components. A prominent example 
would be the problem of spin liquids where one often has a 
description involving matter e.g. spinons coupled to gauge 
fields. Both the matter and gauge field physics can be inter- 
esting, but only together do they form the unified spin liquid 
state. In this context, the prototypical entanglement sum rule 



is the statement that S = S„ 



S g where S„ 



and S come 



from matter and gauge fields respectively. This sum rule is 
a non-perturbative statement about entanglement entropy that 
greatly extends our ability to compute entropies and can be 
very useful for comparing with numerical computations. Such 
a rule has previously been obtained in the context of various 
topological phases (71 Il2ll in a certain limit, but our result is 
much more general. Of course, not every spin liquid state 
admits an entanglement sum rule, but many interesting states 
do. We also emphasize that the notion of an entanglement sum 
rule is general and is not restricted to spin liquids. 

In this work we prove a general entanglement sum rule for 
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a wide class of exactly solvable models that include in var- 
ious limits free fermions, the toric code, and free fermions 
coupled to the toric code. A general model encompassing all 
these limits was recently introduced in Ref. [29] to describe 
so-called orthogonal metals. In the context of the exactly solv- 
able model our sum rule applies to the full entanglement en- 
tropy, but more generally, our sum rule indicates that the uni- 
versal terms in a variety of phases and phase transitions will 
be additive in the sense described above. Note that what terms 
are considered universal depend on the nature of the phase. An 
important advance over the results in Refs. |0, EH] is that we 
can also vary the gauge field dynamics (both previous compu- 
tations were down in an extreme deconfined limit where the 
gauge field does not fluctuate). Below we describe an exactly 
solvable model, introduced in Ref. [29] in the context of non- 
Fermi liquid metals, and prove the entanglement sum rule for 
this model. We also generalize the model while preserving the 
sum rule. Ultimately our sum rule applies to a huge variety of 
physical states including Fermi liquids, non-Fermi liquid met- 
als, spin liquids, deconfined critical points, and much more. 



II. EXACTLY SOLVABLE MODEL 

Consider a square lattice with fermions (electrons) c r on 
sites and Ising spins cr* r , on links. Following Ref. 12911 the 
Hamiltonian is taken to be 

H = —w c^ r G z rr ,c r ' — c\c r 

<rr'> r 

- 9 J ]T J^(- 1 ) 4C " II ^ 

<rr'> r <r'r> 

p <rr'>£p 

The J term is a product over all links sharing a vertex r while 
the U term is a product over all links in given plaquette p. We 
will see that this Hamiltonian describes a (highly fine tuned) 
transition from a Fermi liquid to an orthogonal metal as a 
function of g. Indeed, this Hamiltonian is actually exactly 
solvable for all values of the parameters, but first let us get a 
sense of the physics as a function of the coupling g. 

When g ^ 1 the spins want to polarize as a z rr , = 1. The 
fermions then decouple from the spins and are simply de- 
scribed by a free Fermi gas with hopping w and chemical 
potential fi. In the opposite limit, when g -C 1, the spin 
Hamiltonian is that of the toric code or of Z^ gauge theory 
with gapped matter, and the fermions c r couple minimally to 
the gauge field. Thus in this limit we have again a free Fermi 
surface coupled to the tensionless limit of a Z2 topological 
phase. The topological structure tells us that there has to be 
transition between the g 3> 1 and j<1 phases. 

The U term commutes with everything else in the Hamilto- 
nian and hence we may consider only the subspace in which 

*p = n<rr'>e P a rr' = 1 for a ^ P- In our interpretation 
above, $ p is the Z 2 flux through the plaquette p. Within this 
subspace the a variables are constrained, so we can introduce 



new variables via 

<7 Z r i = T^T^i (2.2) 

and 

IJ a x rr ,=T z r . (2.3) 

<r'r> 

The commutation relations of these operators are preserved by 
this identification provided the r variables also obey the stan- 
dard Pauli commutation relations. Physically, r^f measures the 
Z2 charge of a site while t^t*, creates Z2 charges at r and r' . 
Note that ordinarily we would need to stretch a Wilson line 
between two such charges, but because the flux $ p is exactly 
one, no Wilson line is necessary. 
If we now define 

fr=T?Cr, (2.4) 

r x r = if, (2.5) 

and 

f z r = {-lf lc "r z r = (-1)^7?, (2.6) 

then the f and / variables all commute and the Hamiltonian 

is 

<rr'> r 

-gj J J2 f r- W 

<rr'> r 

Thus we have decoupled fermion (/) and tranverse field Ising 
(f ) systems and hence the full Hamiltonian is solved in terms 
of these two models. We have also dropped the U term since 
it is simply an additive constant within the constrained Hilbert 
space. 

We see immediately that // correlations are identical to 
that of a free fermion model for all g. In particular, since 
n r = c) r c r = f},f r it follows that the physical density-density 
correlator is given by the free fermion result for all g. Thus 
the number fluctuations in a region R, defined as the variance 
of the operator Nr = ~}2 reR n r , is given by the free fermion 
value for all g. Furthermore, the thermodynamics of the model 
exactly factorizes so that 

S(T) = S f (T) + S Ising (T) (2.8) 

where S(T) is the thermal entropy. As a technical subtlety, 
this entropy result only applies if we first send U — > 00 so 
that no states with $ p ^ 1 are excited since otherwise the r 
variables are insufficient to describe the full Hilbert space. In- 
deed, since the Ising part at low temperatures is either gapped 
or at worst a 2d CFT with Si S i ng ~ T 2 it immediately follows 
that the thermal entropy is dominated by the / Fermi surface 
at low temperatures. 
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in. ENTANGLEMENT ENTROPY 

We now turn to the entanglement entropy. How we com- 
pute the entropy depends on what cut we choose to make. For 
example, if we regard c and a as the local degrees of freedom, 
then we cannot make a cut in the t variables since these are 
non-locally related to the a variables. To begin, let us deter- 
mine the ground state |0) of the system. Clearly we have 



|o) = |/K M ))|f( 5 )) 



(3.1) 



where |/) and |f } are the ground states of the decoupled sec- 
tors. To compute the entropy due to a c, a cut we must rewrite 
this state in terms of the appropriate variables. 

As a first step, let us consider the relation between r and f . 
Since f z = (— l) c c r z , we can consider the unitary 



U = (1 - c f c) ® 1 + Jc ® t x (3.2) 



which satisfies 



£/ f = U, (3.3) 



U 2 = 1, (3.4) 



and 



U-r*rf = Ut z U 



= (1 - C f c) (g)T Z + <g> {-T Z ) 

= (-i) c V. 

Using r U r we can convert all r variables into f variables, 
and hence it follows that 




c{w,p))\r{g)), 



(3.5) 



or in other words, the state in the c, r variables differs from 
that of the /, f variables by a local unitary transformation 
(product of CNOT gates). Indeed, we see that if we made a 
cut in the c, r variables then the entanglement entropy would 
obey S = S c + S T since r U r is explicitly a product of 
single site unitary operators. 

The next step is to relate the r variables to the a variables. 
The state \r(g)) corresponds to some particular state \a) as a 
function of the a variables, but we must also account for the 
unitary V = r U r . If we require the number of fermions to 
be even, then every term in the action of V on \c)\a) amounts 
to a product of N c t x or N c /2 pairs of t x operators (N c is the 
total number of fermions). Furthermore, every pair of t x s can 
be written as a string of a z operators and hence every term in 
the action of V on the state is the product of a large number 
of a z operators. Note also that the choice of this product is 
largely arbitrary since the ground state obeys $ p = 1 for all 
p. The density matrix is therefore 



Now consider an arbitrary term in the expansion of V of the 
form 

]J n r T r II ( 1_rl r) ] [ ™r< {} (1 - U r ) 

reC R reR-Cn r£C fi re R-C R 

_(3.7) 

where Cr and denote collections of sites in R and R. As 
above, we can replace products of t x with products of a z , and 
in the case when \Cr\ is even, we can choose the a z config- 
uration to lie entirely within R and R. Similarly, in the case 
when \Cr \ is odd we must have one a z on the boundary con- 
necting R and R. Write |c) = \c, 0) + |c, 1) where 0,1 = N R 
mod 2 labels the parity of Nr (which is also the parity of 
iV^). Within a given parity subsector, we can completely fac- 
torize the action of V (by choosing a suitable pattern of cr z s) 
into unitaries acting on R and R separately, 

V\c,x) = V Rx V Rx \c,x). (3.8) 

This fact enables us to write the state as 

V RO V RO \c,0)\a) + V m V R1 \c,l)\a) (3.9) 

where V x preserves the parity of \c,x). 

Let us now compare the reduced density matrix of R with 
and without the action of V. Without V we have 

tr R ((|c, 0) + \c, l»«c, 0| + (c, 1|) ® \a)(a\) (3.10) 
which reduces to 

tr« ((|0>(0| + |1)(1|) (8, \a)(a\) = + (3.11) 
since the cross terms cancel. Now with V we have 

tr« (Vb|0><0| (8) \a)(a\V^ + V R1 \1)(1\ «8 \a)(a\V R1 ) 

= VroP ( rIv1 + V riP ^V r1 . (3.12) 

Note that once again the cross terms cancel since, for example, 
we can perform the trace in the number basis with no matrix 
elements between different parity sectors. However, we now 
see that since the state of R before V acted was diagonal in 
parity, it follows that the state after V acted can be related to 
the state before V acted by a unitary transformation. Define 
V R to be 



V R0 
Via 



(3.13) 



where the two components refer to the parity sector. Comput- 
ing 



we obtain 



since Pr ^ x has a definite parity of x. To summarize, we have 



^opMo + ^i/4°Mi- 



(3.14) 



(3.15) 



p R = tr c cr ^(|0)(0|) = tr^,(y|c)|cr)(c|(CT|l /t ). (3.6) shown that the action of V on the state \c)\a) reduces, after the 
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trace over R, to a unitary transformation of the density matrix 
of R. Hence all spectral data is unmodified by the application 
of V despite the non-local transformation between a and r. 

Our final result is then S n (R) = S n (R; c) + S n (R; a) since 
the state \c)\t) is factorized. Note, however, that the trace is 
still over the a variables (not the r variables) and hence we 
must write the \r) state in terms of a variables. This transfor- 
mation is critical since in the r variables the state is a product 
at both large and small g, but this is not true for the physical 
a variables. In words, we find that the total Renyi entropy is 
simply the sum of the Renyi entropy of fermions with band- 
structure specified by iu, fi (with no coupling to a) and the 
Renyi entropy of the a model (with no coupling to c). We 
emphasize again that while this means the fermionic compo- 
nent of the entropy may be computed directly from the state 
|/), we must still convert from the non-local t variables to the 
local a variables to compute the a contribution to the entropy. 

The Ising part of the entanglement entropy never scales 
faster than L (L is the linear size of R), and thus the entangle- 
ment entropy has an L In (L) term determined by the / Fermi 
surface for all g. The entropy of the system at g — oo is sim- 
ply given by the / Fermi surface while the entropy at g = is 
given by S = Sf + (\dR\ — 1) In 2 in accord with the results 
of Ref. 10, [l2tl (as usual, we do the splitting over a variables 
as in Ref. [10]). Our argument works at the level of a uni- 
tary transformation and since the parity structure we used in 
our proof is also a property of excited states, we see that the 
full entanglement-thermal crossover function S n (R, T) is also 
exactly additive [26]. However, we emphasize again that we 
must first send U —> oo for the thermal results to hold in this 
simple form. 

Within the Fermi liquid phase we have thus partially con- 
firmed the universality claimed in Refs. H^[ since we have 
exhibited a Fermi liquid state with with variable quasi-particle 
residue in which the entropy is exactly the free fermion result 
as regards the L In (L) term. Furthermore, because the a vari- 
ables are gapped within the Fermi liquid phase, all of the mu- 
tual information calculations of Ref. 113 ill also go through in 
the Fermi liquid state (the generalization of our results above 
to multiple regions is trivial if the a variables are gapped). Re- 
cently, an interesting numerical calculation 113211 of the Renyi 
entropy in Fermi liquids has also largely confirmed the predic- 
tions of JH Hol concerning the universality of the Widom 
formula with small deviations observed only at very strong in- 
teractions. The precise origin of these discrepancies is not yet 
understood. 



IV. GENERALIZED SOLVABLE MODELS 

We may modify the Hamiltonian in Eq. 12. 1 l in many ways to 
produce other exactly solvable models. The simplest modifi- 
cation is to replace the simple nearest neighbor hopping model 
for the c fermions with a more complex bandstructure includ- 
ing terms like <4<v with r and r' two or more links apart. 
To keep the model exactly solvable each such term should be 
augmented with a factor of Il<rr'>G7 a rr' w ^ ere 7 ls an y 
path that begins at r and ends at r'. We can easily check that 



because <J> p = 1 in the ground state the choice of path 7 is 
irrelevant, and using the change of variables 

a z rr , = t*t*, (4.1) 

we see that 

n<ri=« (4.2) 

and hence each c fermion can still be naturally combined with 
a t x to produce an / fermion. Thus we can produce any shape 
Fermi surface we want as well as Dirac cones and other struc- 
tures. 

As an example of the entanglement sum rule in such a 
model, consider the case where the fermions form Dirac 
cones. In this case, the entropy of the fermion component 
is expected to scale as S ~ aL — b with b universal (see Ref. 
[8] for a calculation of b for a disk). Similarly, for j < 1 the 
entropy of the gauge field component is S ~ a 1 L — b 1 with 
b' = In (2). Our results imply a universal term for the coupled 
system of &Dkac + In (2) which should be compared with the 
results in Refs. 10, [H. 

We can also include density-density interactions for the c 
fermions since we have already seen that c^c = f*f. The / 
fermion model may not be so easily solved, but the decom- 
position of the entanglement and thermal entropies still hold. 
Furthermore, we can substitute bosons, call them b, for the 
c fermions without changing the story, and we can of course 
include boson-boson interactions as well. It is even possible 
to add pairing terms for the fermions, again with the appro- 
priate factors of a z . There is one subtlety, however, since if 
a single boson condenses then the Z2 gauge structure is im- 
mediately lost. In this case, our additivity result no longer 
applies. Of course, a fermion (or boson) pair condensate is 
perfectly consistent with the Z2 gauge structure since neither 
carry Z2 charge. 

As an example of the failure of additivity, suppose the bo- 
son wavefunction before the application of V (or really its 
analog in the bosonic case) is simply \b) = (££) r \(3) r where 
b\(3) = f3\f3). Suppose also that the gauge field sector is in the 
g 3> 1 limit with the state \<j) given by a sum over all closed 
loops in the a x basis. The transformation which attaches the 
ends of strings to the bosons (the analog of V) is 

H = ( ^(i±(zl)^ + irf: T .). (4 .3, 

Following the discussion above, if b were not condensed then 
we would divide the wavefunction into even and odd parity 
sectors and the argument would proceed. However, since par- 
ity is no longer well defined due to the condensate, it fol- 
lows that the state of R before and after V& acted differ in 
entropy. This is because, as shown in Eq. 13.121 the applica- 
tion of V partially decoheres the state of R into parity sectors 
thus adding entropy. Indeed, assuming that the boson num- 
ber fluctuates a lot, each parity sector will have roughly the 
same weight and V& will add In 2 entropy (associated with the 
parity bit), exactly canceling the topological contribution of 
— In 2 and indicating that the state is topologically trivial. 
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We can also trivially generalize the model to higher dimen- 
sions. Returning the Eq. 12.11 nearly every term generalizes 
immediately to higher dimensions. Let us take a cubic lattice 
in d — 3 as an example. The fermion hopping terms obviously 
generalize as does the gj term. The appropriate generaliza- 
tion of the J term is to associate with each site r the operator 



n - 

<r'r> 



(4.4) 



where the product is over the six links emanating from r. Fi- 
nally, the U term is unchanged in form, but we must now sum 
over all faces of the cubic lattice. With these simple modifi- 
cations the algebraic structure of Eq. 12. H is preserved and all 
subsequent changes of variable go through as in d = 2. 

Another direction for generalization is to consider a Z n 
gauge field. We introduce a link variable Z satisfying Z n = 
1 and a conjugate variable E which satisfies e 2wlB / n Z = 
e 2Trt/n£ e 2TiiE/n^ % ma y ^ e ^ e unc j ers tood as incrementing 
E by one, while E may be viewed as measuring the electric 
flux along a link. Now consider a square lattice in d = 2 with 
Hamiltonian 



H = —W \_. c\,Z rr 'C r i — /LtVJcJfV 



<rr'> 



-gj J2 (Z rr , + Z; r ,) ~ jJ2e 27rlctcr/n I] e 2 ™^-'/" 



<rr'> 



<r'r> 



-u 



n 



p \<rr'>£p 



n 



zt 



(4.5) 



<rr'>£p 



At large g we have Z n gauge theory in the tensionless limit 
or the equivalent string net model while at small g the model 
favors Z rr i = 1 on all links. Repeating our analysis above 
with parity replaced by the fermion number mod n gives an 
essentially identical entanglement sum rule. There is, how- 
ever, a greater constraint on the allowed fermion states, since 
we must leave the Z n gauge structure unbroken. The sum 
rule will break down for fermion pair condensates but will re- 
main valid provided we condense a multiple of n fermions (or 
bosons). 



V. DISCUSSION 

In this work we have derived an entropy sum rule for a wide 
class of exactly solvable models. Furthermore, it is clear that 
this technique has not been exhausted and may give further 
useful insights into the entanglement structure of quantum 



matter. There are many simply extensions of our arguments 
here, some of which we have already sketched, but it is proba- 
bly also true that other classes of models admit structures like 
the entanglement sum rule. 

For the exactly solvable models we considered, we have 
proven exact additivity for the full entropy. More generally, 
however, our results imply that the universal terms will add 
in any phase smoothly connected to one describable by our 
exactly solvable models. These phases include the Fermi liq- 
uid, the orthogonal metal 12911 . various kinds of spin liquids, 
exotic superconductors, orthogonal nodal metals, and much 
more. Moreover, we have the added ability to follow the en- 
tropy through various quantum critical points, and provided 
the critical point is not heavily fine tuned (as unfortunately 
is the case with the orthogonal metal), universal terms at the 
critical point should add as well. 

We also give an example of a broader application of our 
ideas. It was noted in Ref. [29] that Luttinger's theorem, when 
take to refer to singularities in the electron spectral function, 
is not satisfied in the orthogonal metal e.g. the electron spec- 
tral function sees an overabundance of critical surfaces rela- 
tive to the / Fermi surface. Of course, Luttinger's theorem 
only applies in its original form to Fermi liquids, so nothing is 
surprising about the observation of Ref. 12911 , but it has been 
possible to preserve the counting embodied in Luttinger's the- 
orem in other phases. Here we see that the Luttinger count 
may be upheld provided we use the entanglement entropy to 
measure the hidden / Fermi surface. This is intriguing in that 
it makes contact with recent attempts to measure hidden Fermi 
surfaces in compressible states using entanglement entropy in 
a holographic context 13411 . What we have shown is that the 
entropy does exactly capture the hidden Fermi surface in the 
compressible orthogonal metal. 

Finally, let us conclude by noting that many recent experi- 
ment measurements on the organics l3~5ll36ll and other two di- 
mensional magnetic systems have revealed evidence of a spin 
liquid ground state, possibly one including a Fermi surface 
of spinons. However, the identification of the precise spin liq- 
uid state has remained elusive with various experiments point- 
ing in different directions. Thus more detailed information is 
needed, and it may be that the detailed information on entan- 
glement presented here when combined with numerical work 
can help elucidate the precise nature of these exciting new 
phases. 
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